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1. Introduction 

The study of spin systems on finitely connected random graphs started nearly thirty years 
ago [U El El HI [5] , but has in the last decade enjoyed renewed popularity as a result of many 
successful multi-disciplinary applications of the mathematical tools that it generates. The 
reason is clear: interacting particle systems in the real world do not have full connectivity, 
but generally involve an average number of interaction partners per unit that is indeed 
limited. Moreover, apart from carefully prepared pure samples of magnetic materials or 
crystalline solids, the graph that represents which elements interact with each other appears 
in most disciplines random in first approximation, often with features that are remarkably 
universal across fields, such as a power law distribution of node degrees. The mathematical 
techniques used to study such systems are consequently being refined and applied in areas 
as diverse as spin-glass and glass modeling [SI [TJ El El QUI Uj, error correcting codes 
[El 00 EE [T6], theoretical computer science pH (H QH [201 ED El [23] , recurrent neural 
networks [M1E3ES], 'small- world' networks [271 ESI EFJ, socio-economic modeling [301 ETJ, an d 
CDMA communication [321 E3J [Ml 155] . Following a first wave of equilibrium studies, based 
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on the finite connectivity replica method, we are now also beginning to acquire mathematical 
techniques with which to bring the dynamics of spin systems on finitely connected random 
graphs under control [36l E3 EH ESI SHI [HI S21 H3]. 

In modeling real-world systems of the above type one places the dynamical variables 
(e.g. spins) on the vertices of a graph, with edges connecting the variables that interact. 
The behaviour of the model depends strongly on the choice made for the graph; however, 
this graph has to be drawn randomly to make analytical progress using disordered systems 
theory. So one tries to define a random graph ensemble with topological characteristics that 
are as close as possible to those observed in the system to be modeled, while keeping the 
model solvable. The state of the art in this respect is defining a random graph ensemble where 
edges are drawn subject to the constraint that all vertex degrees are prescribed. Here the 
only information on the topology of the system that is effectively carried over from the real 
world (apart from irrelevant site permutations) is the graph's degree distribution. Since for 
most degree distributions there is still a large and diverse set of compatible microscopic graph 
realizations, with possibly distinct macroscopic phenomenology for spin dynamics and statics, 
one would like to increase the amount of topological information embedded in the random 
graph ensemble, beyond prescribing just the degree distribution. 

In this paper we present a simple class of Ising spin models on finitely connected random 
graphs, where these graphs are drawn from ensembles in which not only the degree distributions 
can be chosen arbitrarily and imposed as a constraint, but where in addition the edges are 
drawn in a way that allows for preferential attachment on the basis of an arbitrary function 
of the degrees of the two vertices concerned (similar in spirit to the so-called hidden variable 
ensembles |444 l4"5j l4"6l W7\ 148] ). The graphs thus generated are no longer characterized by 
their degree distribution alone, yet the associated spin models can still be solved in thermal 
equilibrium using conventional finite connectivity replica techniques. The graphs in the 
proposed ensemble remain effectively tree-like, and the clustering coefficients of randomly 
selected vertices are zero with probability one (as in the standard finite connectivity ensembles 
with prescribed degrees only). We solve our models with the replica symmetry (RS) ansatz, 
and show which features of the phase diagrams and the dependencies of observables on control 
parameters are identical across all graphs with the same degree distribution, and which features 
depend more specifically on the extra topological information that is included. 

2. Definitions 

We study a finitely connected and bond-disordered system of N interacting Ising spins 
crj £ {—1,1}, in thermal equilibrium characterized by the following Hamiltonian 

H ( a ) = ''■j (J '- , 'J C7 i cr = (di, . . . , a N ) (1) 

i<j 

The frozen variables Cj,- £ {0, 1} define the connectivity of the system; they define a random 
graph, with vertices labeled by i,j £ {1, . . . ,N} and with Cj,- = 1 if an only if i and j are 
connected by a link. We define Cy = Cji and cu = for all (i,j), and abbreviate c = {cy}. 
The bonds Jij £ IR are drawn randomly and independently from some distribution P(J). To 
characterize the topology of a graph c we define for each vertex i the degree h L {c) = ■ c ij (the 
number of links to this vertex), and the degree distribution p(k\c) = N^ 1 J2i <5fc,fc;(c)- Thus the 
average connectivity of a graph c is (k) = J2k>o kp(k\c), which we choose to be finite, even in 
the limit N — > oo. We draw the graph c randomly from an ensemble defined by a probability 
distribution Prob(c) in which not only the degrees are constrained to take prescribed values 
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{ki, . . . , kjy}, but where the link probabilities are modified further according to some function 
Q(., .) of the degrees of the two vertices involved: 

Prob(c) = ^~li [^Q(ki,fci)<W + (l—Qfak^J II<W) (2) 



N J 3 V J\f 

- i<j - s 

z * = EI1 [^Q(^,%)<W + (i- { -^Q(k,k 3 ))5 Ci J ]j8 kiMc) (3) 

The A?" degrees are, in turn, drawn randomly from a prescribed distribution p(k). Clearly we 
require Q(k, k') > for all k, k! . In order to ensure furthermore that for large N such graphs 
can actually be found, we need to choose the function Q(.,.) such that in formula ^) the 
partial measure Ili<j [• • •] is consistent with the average connectivity (k) = J2k kp{k) imposed 
by the constraining factor FJj ^^(c)- Upon defining Q(k, k') = Q(k', k), this is achieved when 
lirn/v^oo iV -2 Q(ki, kj) = 1, so the function Q(., .) is to be chosen subject to 

Q{k, k') > VA;, k' and ^ p(k)p(k')Q(k, k') = 1 (4) 

k,k'>0 

We note that an alternative but mathematically equivalent way to write the graph probabilities 
is obtained by applying to ([2]) the general identity 

IL-><A,.: • = (n^)e E -^- [losA ^ losBd (5) 

which gives 



Prob(c) = — QE^^jU^mQih^/m-logli-i^Qik^/N]} -q 5k ^ ^ 
Z N . 

Zm ^^T, 1<j ^A lo d(k)Q(h,k J )/N]-log[l~{k)Q(k l ,k J )/N]}YT 6 (7) 

c i 

Given the above definitions, our objective is to calculate for the system ([TJ the asymptotic 
disorder-averaged free energy per spin /, in order to find the phase diagrams for spin systems 
defined on typical graphs in the ensemble ([2]): 

1 : — — ^ttt^t , , 1 



/ = - lim — logV e-Wlv) = - lim lim log [V e-^^V (8) 

J N^oo (3N ^ Ar^oon^O f3nN i w 

in which 7TT denotes averaging over the disorder in the problem, viz. the randomly drawn 
graph c with statistics ([2]) and the random bonds {Jij}- This calculation is done with the 
finite connectivity replica method, within the replica symmetric (RS) ansatz. We will be 
particularly interested in the dependence of the phase diagrams on the choice made for Q(., .). 
In the absence of the degree constraints, the function Q(., .) would have controlled the bond 
probabilities fully, via (cy) = Q(ki,kj)(k)/N. Here, in contrast, its role is to deform the 
measure imposed by the degree constraints, biasing the probabilities in those cases where 
there exist multiple graphs with the same degree distribution. 

A relevant question to be asked at the start is whether the ensemble deformation induced 
by Q(k, k') can have sufficient impact in the thermodynamic limit on macroscopic observables 
to justify the present calculation. For instance, upon calculating for the ensemble ([2]) the 
joint distribution of degrees and clustering coefficients for — > oo one finds that, like 
ensembles with degree constraints only, the clustering coefficients are zero with probability 
one, see Appendix A However, upon reflection it becomes clear that the proposed ensemble 



deformation will generally affect the system's phase diagram. A quick way to see this is to 
compare the two choices Q(k,k') = 1 and Q(k,k') = 5kk' I J2k" P 2 (k")- In the first case we 
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return to the degree-constrained ensemble in In the second case our ensemble describes 
graphs that are each composed of disconnected regular sub-graphs, of sizes p(k)N for each 
k with p(k) > 0. The transitions away from the paramagnetic state will now be those 
corresponding to the sub-graph with the largest degree allowed by p(k); these will only coincide 
with those of Q(k, k') = 1 when p{k) = 5k, c , i-e. when the degree distribution itself is that of 
a regular graph. 



3. Equilibrium replica analysis 

3.1. Derivation of saddle-point equations 

As usual, we calculate ([8]) upon writing the Kronecker 5s of the degree constraints in integral 
form, using 5 nm = (2ir)~ l J^ n duj e luJ ( n - m ) , This gives, after some rearranging of summations, 
factorization over the disorder variables. We define the short-hands <jj = (aj, . . . , <r"), so that 

/= Jim lim. — \\ogZ N -log / TT l"^e" 



N^oo rwO (3nN 



CT i...(T jy I 

II (l + JfQih, kj)[J dJ p(j) e ^-^-i(^+^)- l]) } 

i<j 

Jim Mm _L{ \ og z N - log Y r TT [^ e iw * 

(J i ...CT jy ) 



V 



x exp 



l -§fJLQ( k ^ k j)l[ dJ J P(J)e /3J<Tl ^- i( ^ + ^ ) - 1] +O(iV )]} (9) 

y 

We proceed towards a steepest descent integration by introducing for <x £ {— 1, l} n and 
A;G {0, 1, 2, ...} the functions cr|{<Tj, Wj}) = A r_1 5k,k i 5tT,cr i e~ luJi - They are introduced 
via the substitution of integrals over appropriate ^-distributions, written in integral form, viz. 



dD(k, cr)dD{k, a) i ND(k,cr)\D(k,cr)-D(k,cr\{cr i ,w i })] 
2-k/N 



(10) 



Upon using also iV -1 J2ij Q{k, %) = N + O(VN) due to (gj), and the short hand {d£>dZ?} = 
<j D(k, cr)dD(k, tr) we then obtain 

/= l im lim-L/log^-log /{dI)dD}e i7V S^^' T ) D ( fc '' T )-^W +0 ^ 1/2 ) 

xexp [-(k)NyQ(k,k') y D{k,a)D{k',a') fdJ P{J)e 
kk' acr> 

x exp \Ny v (k)\o^y [* ^-w^ i 



J3J<T-(T' 



We next define z = limjv^oo iV 1 log-Z/v (anticipating this limit to exist), which allows us to 
evaluate / by steepest descent and write 

1 

— i 

n 



f 



= l^ - eXtl {D,D}fn[{ D ^ D }} 

f n [{D,D}\ = - 1 ji^^ffW,^) - \{k)~z 



(12) 
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+ £ p(k) log £ f ^e-*-^'^-" ) (13) 

The extremization in (|12p with respect to {-D,-D} gives the following saddle-point equations: 
D(k, cr) = i(Jfe) ^ Q(Jfc, fc') J] cr') / d J P( J)e /3Jcr - cr ' (14) 

fc' <X' 

p(fc) r do; e ^( fc - 1 )- i ^( fc ' cr ) e " i " 
a) = V ;J ~^ s — (15) 

The second of these equations can be simplified using the identity 

r du; e ^D(k,a)e-^ = f 2tt[-LD(A;, if I > 

»/ — 77 \ if £ < 

So, if we also re-define D(k,cr) = i{k)F(k,cr), we arrive at 

F(k, (t) = Y, Q( k , k ') S W / dJ p (J)^ J<7 ' a ' (17) 
fe' cr' ^ 

p(fc)fc F k -\k,a) 

D{Ka) - ^T^Hk^) (18) 

3. ,2. Simplified expression for free energy per spin 

Formula (|12|) for the disorder-averaged free energy per spin, which is to be evaluated at 
the relevant solution of the saddle-point equations (|17ll8p . still contains the term z = 
liniTv— >oo N ~ 1 log Zjy, which measures the effective number of graphs in our ensemble ([2]). 
Since z is independent of f3 we can use the identity lim v g_ > o(/3/) = — log 2 to find it. With (fl"6j) 
and (|17I18P we first write (HU as 

7 = }^j-{z + (A) - (k) log(fc) - log [1 £ } (19) 

(with D and F taken at the relevant saddle-point). Working out the saddle-point equations 
for (3 — > shows that there the two order parameter functions {D, F} become independent of 
cr, viz. D(k,cr) = 2~ n D{k) and F(k,cr) = F(k), where the latter obey 

F (k) = E P ~^Q^ ^(k'), D(k) = ^F-\k) (20) 

Insertion into the equation lim j g^o(/3/) = — log 2, together with (|19p . then leads us after some 
further simple manipulations to the following formula, where F{k) is the solution of (j20[) : 

z = (k) log(fc) - (fc) + X; P(k) log (21) 

and hence 

7= " lim -^5:p(A ; )log[5:[F(A ; , < T)/F(A ; )] fc ] (22) 
' fe cr 

We note that for the non-deformed graph ensemble with degree constraints and finite 
connectivity statistics only, viz. Q(k,k') = 1 for all (k,k ! ), the solution of ([20]) would be 
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F(k) = 1 and D(k) = p(k)k/ (k) . If we now define n(k) = e ^(k) k /k\, i.e. Poissonnian 
degree probabilities with average degree (k), we see that 

z =z nd + Y,P(k)k log F(k) (23) 

k 

Znd = E K*0 l0 § < k ) = ~ H P ~ D (P\ K) ( 24 ) 
fc 

with the entropy H p = — J2kP(k) log p(k) > of the degree distribution and the Kullback- 
Leibler distance D(p\\ir) = J2kP(k)log\p(k)/ir(k)] > between the actual degree distribution 
p(k) and the Poissonnian ir(k). Since z measures the effective number of graphs that can be 
generated from our ensemble, and z n d is its value in the absence of deformation, it follows 
from (|23p that we can define a simple measure Ed e f of the graph specificity increase resulting 
from the introduction of the deformation defined by a function Q(k, k') as follows 

Edef = -5>(A:)fc log F{k) (25) 
k 

3.3. Replica symmetric theory 

To take the required limit n — > in our formulae we make the ergodic or replica-symmetric 
(RS) ansatz. The replica order parameter D(k,cr) must now be invariant under all replica 
permutations, and therefore have the following form, with J dh D(k, h) = J2a F>{k, 

We work out the implication of this ansatz for the order parameter F(k, a), using the identity 
f(cr ) = e Aa B, with A = \ log[/(l)//(-l)] and B = V/(l)/(-l) (which holds for a = ±1), as 
well as the identity | log[cosh(x + y)j cosh(x — y)\ = atanh[tanh(x) tanh(y)]. This results in 

F(k, *) = y: q% I h')p(j) n cosh ^«+ h 'v 



with 



J2Q(k,k') dh'dJ ^(fc>0^(^)Gn(/iU)e5 ( S^" )los[cosh ^ [J+ ^ )/cosh(/3[J - h ' 1)] 
k' J 

JdhF(k,h) (27) 

G n (h,J) = [cosh( ( 5[/i+J])cosh(/3[/i-J])/cosh 2 (/3/i)] n/2 (28) 

F(k, h) = Y,Q(k, k') f dh'dJ D(k', ti)P(J)G n (ti, J) 
k' J 

x 6[h - /? _1 atanh[tanh(/3J) tanh(/%')]] (29) 



We will now have two new saddle-point equations, written in terms of the RS kernels F and 
D. The first equation is (|29p . The second follows upon inserting (|26l27p into (| 1 S |) : 

dh D{k h) J h Z a ** _ P (k)kSIU<k-ShtF(k,h i )]^ E ' <raEl ^- lh ' 



[2cosh(/%)] n (k) JU e <k[dh,F(k,h e )][2co S h((3j:i<khi)] r 
p(k)k 
= lkj 

JT\e<k~i[dheF(k, h e )]5[h - J2e<k-i he] 



. Jdh fP h T, a °« 



IUe<k [dh e F(k, h e )] [2 cosh(/3 J2 e<k h t )]' 



(30) 
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From this result we can read off the second order parameter equation. Upon taking the limit 
n — > in the latter result and our first equation (|29p . we arrive at the transparent expressions 

F(k, h) = Y^Q(k, k') f dh'dJ D(k\ h')P(J)5[h - -atanh[tanhQ3J) tanh(/3/i')]] (31) 
k' ■* " 



D(k,h) 



P (k)k IUe< k -i[dheF(k,h e )]5[h-E 



Kk-\ 



he 



(32) 



(k) [Jdh'F(k,h')] k 

Upon defining finally D(k) = Jdh D(k, h) and F(k) = JdhF(k, h), we discover that these last 
two integrals are exactly the quantities in (|20p . since upon integrating over h in both equations 



(|31|32p they reduce precisely to (1201) . This does not seem to be a trivial result, since (j20|) was 
derived from the properties of the random graph ensemble alone. Given these relations one 
is prompted automatically to define D(k,h) = D{h\k)D(k) and F(k,h) = F(h\k)F(k), where 
now J dh D{h\k) = Jdh F(h\k) = 1. Our RS order parameter equations thereby take the form 

F{h\k) = Y iS^T^mTK fdh'dJ D{h'\^)P(J)8[h~a.t^[tBjik(fiJ)t^{ph')^ (33) 
k i {k)F{k)F(k ) J p 



D{h\k)= H[dhtF(h e \k)]6[h-J2h e ] 



£<k 



F(k) = (ky 1 J2 P (k , )k , Q(k,k')F- 1 (k') 



(34) 
(35) 



where D(k) subsequently follows from the second identity in (|20p . and where ()35[) automatically 
ensures that the solutions of (|33|34p are normalized. If we write the free energy (|22p for our 
present RS solution in terms of the quantities in (|33|34|35p and take the limit n — > 0, using 
relations such as J2k D(k)F(k) = 1, we find in a similar manner the remarkably simple result 



/ 



RS 



\ lo § 2 " \ E K*0 I II [dhtF(h t \k)] log cosh(/3 Y, h i) 

P P k J Kk Kk 



(36) 



3.4- Physical observables 



To find formulae for observables like m = Mm^-^oo iV -1 J2i an d Q = lhnjv-voo N^ 1 J2i ( a i) 2 i 
it will be convenient to calculate for an re-replica system where <jj = (aj, . . . , erf), i.e. before the 
limit n — > 0, the expectation value P(k,a) = limN^oo N^ 1 J2i (<Sfc,fc i <S<r ) <r j )- We can calculate 
P(k, <r) using steps similar to those taken in the evaluation of the free energy, using (|13|) : 



P(k,cr) = lim iV^ 1 V 



a jv 



-PE a H ( aa ) 



lim lim y S (a ctHct1 ^e-^S'^) 



-n N^oo N 
1 



cr 1 rr 



n-\-n f 



lim lim -Y^k t {dDdD}e 

'— »— n Af— >oo IS 1 



-PNf n+n ,[{D,D}} 



Ea n+1> ...,a n+nl /da; e 



iuk-iD(k,a 1 ,...,o n+n ,)e lu 



iwk— iD(k,cr'i ,...,(?' ,)c~ iu - 



1 n+n' 

V fduj J^+(k)F{k,a 1 ,...,a n+n ,)e- i ^ 

Z^cr n+1 ,...,<7 n+n , J 



lim p(A;) 



n' — >— n 



E<,... CT ' ,/dwe 

1 T). + 77.' 



iwfc+<fc>F(fe,eri,...,<7' )e-i" 



(37) 
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where in the last line we now have to take F(k, <7i, . . . , cr n+n i) at the saddle-point of / n+n /[. . .], 
with D[. . .] = i(k)F[. . .], and where we have used lim n ^o extr| D ^f n [{D, D}] = 0. Once more 
we can carry out the ^-integrations and find 

F k (k,a 1 , . . . ,<7„ 



P(k, a) = lim p(k) 



y 



>n+n' , 



(38) 



E a , , F^k,a[,...,a' n+nl ) 

1 n+n' 

In replica symmetric states, where F(k,cr) = Jdh F(k)F{h\k)e^ h ^<=' eTo! , we can carry out the 
remaining spin summations, giving 

p (k ) y a, [ /n^< fc [d^F(^|fc)][2cosh(/?^< fc ^)]" , e^^ M(E ^"' TQ) 
RS ^ >*> nl 1 ^ > [ mt< k [^F(h f \k)][2 C osH^< k hd] n+n ' 

f p/ 3 (E«<fe^)(Ea<»°'«) 



p{k) dh W(h\k), 



[2 cosh(/3/i)]« 

with the degree-conditioned effective field distribution 



W(h\k) 



l[[dh e F(h e \k)] S[h-J2hi 

Kk ' Kk 



The order parameters m and q can now be written in their usual form 

dh W(h) tanh 2 (/?/i) 



m = Jdh W{h) tanh(/3/i), q 
W(h) = Y,P{ k W{h\k) 

k 

and the free energy formula (|36|) becomes simply 



/ 



RS 







log 2 



dh W(h) logcosh(/3/i) 



(39) 
(40) 

(41) 
(42) 

(43) 



3. 5. Simple solutions for special cases 

The simplest limit is /3 — ► 0, where we expect a paramagnetic state. Setting /3 = 
in our order parameter equations (|33l34l35p indeed gives F(h\k) = D(h\k) = 5(h), and 
limfl_*o(/3/) = — log 2 as well as m = g = 0. As always in such systems one notes that 
F(h\k) = D(h\k) = 5(h) is in fact always a solution of (|33l34l35p . at any temperature. 

A less trivial special case is the choice Q(k, k') = 1 for all (k, k') (non-deformed graph 
ensemble), where we should be able to connect the present theory to earlier results in literature. 
Here one has F(k) = 1 and D(k) = p(k)k/{k), and F(h\k) is no longer dependent on k. We 
may thus simply write F(h\k) = F(h) (not to be confused with F(k)), so that upon eliminating 
D(h\k) from (|33[) via (|34p . we are left with the order parameter equation 



F(h) 



\ - p(k)k 



dJ P(J) /jJ[d^ J F(^)]j[^--atanh[tanh( ( aJ)tanh(/3^/t £ )] 

J Kk ' Kk 



(44) 



To establish the connection with earlier results for the Ising system on graphs taken from the 
non-deformed ensemble we define a new field distribution W(h) 

p(k)k 



W(h) 



H[dh e F(h e )]5[h-Y l hi 

Kk ' Kk 



(45) 
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According to (|4Up . the latter cavity field distribution W(h) is expressed in terms of the 
degree-conditioned effective field distributions W(h\k) via W(h) = J2k p(k)k(k)^ 1 W(h\k — 1). 
Equation (|44p now tells us that 

F(h)= JdJdh'P(J)W(h')5[h-^a,tanh[tanh{pJ)t&nh(ph')}^ (46) 

and hence we find the following RS order parameter equation in terms of W(h), which we 
recognize from earlier studies on Ising systems with random graph ensembles that are given 
prescribed degree distributions: 

W{h) = £ / II [dJedh e P(J e )W(h e )]5\h - ~ £ atanh[tanh(/5 J t ) tanh(/3^)] 1 (47) 

k W J Kk Kk 

Similarly we can write the free energy (I36p for non-deformed ensembles in terms of W(h): 
7 RS = -~lo g 2-±5>(*0 J]J[dh e dJ e P(J £ )W(h e )} 

k £<k 

x log cosh ( ^2 atanh[tanh(/3 J e ) tanh(/3/i^)]) (48) 



t<k 



We should emphasize that the two field distributions W(h) = J2kP(k)W(h\k) and W{h) = 

fc) are generally different. The obvious exceptions are Poissonnian 
degree distributions, where p(k + l)(k + 1)(A;)~ 1 = p(k), and systems where W(h\k) does not 
vary with the degree k. 



4. Continuous phase transitions 

4-1. Bifurcations away from the paramagnetic state 

Continuous bifurcations away from the paramagnetic (P) state F(h\k) = 5(h) are found 
in the usual manner, by expansion in moments of F(h\k). We assume the existence of a 
small parameter e with < |e| <C 1 such that fdh h l F(h\k) = 0(e ). Let us first define 
€k = fdh hF(h\k). Multiplication of (|33|34p by h, followed by integration over h gives the 
lowest nontrivial order: 

£ fe = \ E %f^)F(V) I dJ P{J) /ll ^(^l^)]atanh[tanh^j)tanh^^ ^)] 

= JdJ P(J) tanh(/3J) £ ^^F^ ** + °<?) ^ 

Thus a continuous transition to a ferromagnetic (F) state occurs when the matrix with entries 
Mkk' f dJ P(J) tanh(/3J) has en eigenvalue one, where k, k' G {0, 1,2,.. .} and where 

Q(k,k')p(k')k'(k'-1) 
Mkk ' = (k)F(k)F(k>) (50) 

In a ferromagnetic state one will have a nonzero magnetization m = (3J2kP(k)k e k + 0(e 2 ). 
Similarly we can check what happens if fdh hF(h\k) = for all k, so m = 0, and the first 
nontrivial order to bifurcate is e 2 . This corresponds to a transition from a paramagnetic to a 
spin-glass (SG) state. Now we define e k = fdh h 2 F(h\k) = 0(e 2 ). Multiplication of (1331341) 
by h 2 , followed by integration over h now gives 



e k = fdJ P(J) tanh 2 (/3J) £ Q ^^' J \{ [dh e F(h e \k')] £ h e h' e + O 



(k)F(k)F(k' 



e 3 ) 



Kk' ll><k' 
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6J "I^IW£ W w y + Otf (51) 

(where we have used Jdh hF(h\k) = to eliminate terms with i ^ £'). Thus a 
continuous transition to a spin-glass (SG) state occurs when the matrix with entries 
M kk ,fdJ P( J) tanh 2 (/?«/) has an eigenvalue one. The previous P— >F transitions mark 
bifurcations away from (m, q) = (0, 0) to states with (m 7^ 0, q 7^ 0), and the P— >SG transitions 
mark bifurcations away from (m, q) = (0, 0) to states with (m = 0, q 7^ 0): 

P -► F : M kk'£k> = A F e fc Ap 1 = fdJ P(J) tanh(/3J) (52) 

k' J 

P - SG : £ M **' e *' = A SG£fc Asg = f dJ P ( J ) tanh 2 (/3J) (53) 
k' 

with the matrix elements M kk > as given in (|50p and with F(k) to be solved from f)35[) : 

= (ky^Pi^k'Qik, k')F- l {k') (54) 
fc' 

The physical transition occurs at the largest eigenvalue of the matrix M. It will be of the type 
P^F if fdJ P{J) tanh(/3J) > JdJ P(J) tanh 2 (/3J), and otherwise it will be P^SG. If for the 
bond distribution we choose the binary form P(J) = 5(1 + n)5(J — Jo) + ^(1 — n)5{J + Jo), 
we can already deduce that the triple point (where the phases P, F and SG come together) 
is found along the line n = tanh(/3Jo), with the P— >F transition being the physical one for 
n > tanh(/JJo) and P^SG transition being the physical one for n < tanh(/3Jo). This will be 
true irrespective of the choice made for the ensemble deformation function Q(k, k'). 

4-2. Analysis of the eigenvalue problems 

We will now analyze the bifurcation eigenvalue problem for the matrix (|50p . to be solved in 
conjunction with (|54p . for a number of graph ensembles, which in this paper are characterized 
by an ensemble deformation function Q(k, k') and a degree distribution p(k). The deformation 
function must always obey Q(k,k') > 0, Q(k,k') = Q(k',k), and J2kk' P(k)p(k')Q(ki k') = 
1. Upon writing the largest eigenvalue of the matrix ([50]) as A max (Q,p), the continuous 
bifurcations away from the paramagnetic state occur for 

P - F : 1 = A max (Q,p) J dJ P{J) tanh(/3J) (55) 
P -» SG : 1 = A max (Q,p) j dJ P{J) tanh 2 (/3J) (56) 

Type I: Separable deformation functions. The simplest family of deformation functions 
are of the separable form Q(k,k') = g(k)g(k')/(g} 2 , with g(k) > for all k, and (g) = 
J2kP(k)9(k) > 0- The special choice g(k) = 1 gives the ensemble without deformation. For 
this family it follows immediately from (|54p that F(k) = g(k)/{g), and (|50p reduces to 

M kkl =p(k')k'(k'-l)/{k) (57) 

There is just one eigenvector, namely e k = 1 for all k, with eigenvalue A = (k 2 ) / (k) — 1. Hence 
the continuous transition lines are 



F : 1 = [{k 2 )/(k) - 1] / dJ P(J) tanh(/?J) (58) 
SG : 1 = [(k 2 )/(k) -l]fdJ P(J) tanh 2 (/3J) (59) 
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The function g(k) has dropped out of our equations, so the transition lines will be identical 
to those of the non-deformed ensemble, i.e. to those found when Q(k,k') = 1. For a 
Poissonian degree distribution one has (k 2 ) = (k) 2 + (k), so A = (k) and we recover the 
standard results for Erdos-Renyi graphs. Moreover, the solution F(k) = g(k)/(g) of (|54p gives 
Q(k,k')/F(k)F(k') = 1 for all (k,k'), which implies that g(k) also drops out of the order 
parameter equations (|33|34p . It follows that for separable deformation functions Q not only 
the transition lines, but the complete solution of the model, including the values of the physical 
observables anywhere in the phase diagram, is independent of g(k) and therefore identical to 
that of the ensemble with degree constraints only. This is true for any degree distribution, 
and has a simple explanation: upon substituting Q(k,k') = g(k)g(k')/ (g) 2 into the graph 
probabilities © one obtains 



Prob(c) 



iV{i<fc>log[<fc>/<g> 2 7V]+7V- 1 ^ i fc,(c)logg(fe l (c))+0(iV- 1 )} 
3 J [li d ki,ki(c) 



7V{I(fc> log[(fc>/<g) 2 iV]+JV-i fei(c') logs^cOHO^- 1 )} 
L~iq! ^ 11? 

e <"°>n«W) 



EceWn^MC) (60) 
For Q(k,k') = g(k)g(k')/ (g) 2 the Q-dependent factors in Prob(c) depend in leading order on 
c via the degrees fcj(c) only. Since the degrees are constrained, these factors drop out, leaving 
only subdominant terms with a vanishing impact on the thermodynamics in the limit ./V — > oo. 

Type II: Additive deformation functions. The second class of deformation functions we 
will study is Q(k,k') = [g(k) + g(k')]/2{g), with g(k) > for all k and (g) > 0. Again the 
simplest choice g(k) = 1 gives the non-deformed graph ensemble. Now it follows from (|54p 
that F(k) = Ag(k) + B, with 

A _j_ / fc \ p 1 / fcfl(fc) \ ffm 

2{k){g)\Ag(k) + B/ 2(k)(g}\ Ag(k) + B / { ' 

We can rewrite B to get a simple relation between A and B: 

AB = 1 / M \ = _ \ = J_ _ Afl (62) 

2{k){g)\Ag(k) + B/ 2{k){g)\ Ag(k) + B> 2(g) { ) 

Thus AB = 1/4(5), i- e - B = l/4A(g). Upon eliminating B from our equations and upon 
defining A = x/(g) with x > 0, we then find that x is to be solved from J-(x) = 1, where 

HX) = ^(l + Jlw/ig)) (63) 
We note that dJ 7 (x)/dx < for x > 0, with .F(O) = 2 and ^"(00) = 0, so there is indeed 
a unique and well-defined solution x > of J-(x) = 1. For the trivial case g(k) = 1 (i.e. 
Q(k,k r ) = 1, no ensemble deformation) we obtain ^(x) = 2/(1 + Ax 2 ), giving x = 5 and the 
correct simple solution -F(fe) = 1 encountered earlier. 

We proceed with the analysis of nontrivial choices for g(k). Let us define the short-hand 
G(k) = g(k)/(g), so (G(k)) = 1, Q(k,k') = \[G{k) + G(k% and F(k) = xG(k) + l/4x. The 
matrix (|50p to be diagonalized then takes the following form: 

= 8x 2 [G^ + ^fcQMfcQ^^-l) 

(ib) [4x 2 G(k) + l][4x 2 G(k>) + 1] 1 ' 

Its eigenvalue equation becomes (with brackets denoting averages over the degree distribution): 
8x 2 G(k) / k'{k'-l)e kl v 8x 2 / GCfcQfc^fc^-^efc/ v 

fc (fc)[4x 2 G(fc) + l]\4x 2 G(A;') + l/ (fc)[4x 2 G(/t) + l] \ 4x 2 G(/fc') + l ' 
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We see that the components e/% of any eigenvector must always be of the form 

cos(</))G{k) + sin(» 
€k = 4x 2 G(k) + l (66) 

where 

. ^ 8x 2 / k(k-l)[cos(cj ) )G(k)+sm( ( j ) )} \ 

Acos(^) = w ( [4x 2 G{k) + l] 2 ) (67) 

. 8x 2 / G(fc)fc(A:-l)[cos(^)G'(fc)+ S in(0)] v 

AsmW = [4* 2 G(fc) + l] 2 > (68) 

or, in matrix form: 

/ /JL\\ Q 2 / / fc(fc-l)G(fe) v / fc(fc-l) v \ / x 

/ cos(0) \ _ 8s_ [ ^[4i 4 G(t)+ip' I f cos(0) \ , > 

1 sin(0) ) (k) / fc(fc-i)G 2 (fc) y / fc(fc-i)G(fc) v I sin (0) I ^ oy ^ 

The two eigenvalues are now calculated easily. We need the largest one, giving us the following 
expression for X ma , x (Q,p) for the present family of graph ensembles, with y = 4x 2 /(g): 



x fn ^- 2y J/ fc(fc-l)g(fc) \ // fc(fe-l) \ / fc(fe-l)ff 2 (fc) \\ 
maxW ' pj " <*) \\ [ y5 (fc)+i] 2 / + V W)+iF \ [ysW+i] 2 // ( j 

where y is the solution of 

>yg{k) + l> 2 

It is a trivial matter to check that for the simple choice g{k) = 1, where y = 1, one indeed 
recovers from (|70p the correct eigenvalue (k 2 )/(k) — 1 of the non-deformed ensembles. It will 
also be clear from (|7Up that for the present non-separable family of ensemble deformation 
functions Q(k,k'), the phase diagram will generally indeed be affected by the deformation. 

Type III: Simple binary deformation functions. Our third class of deformation functions 
are those where Q(k,k') takes only two values. Here the deformation can even strictly 
forbid links that otherwise would have been allowed. We will focus on the simple example 
Q(k,k') = 7o + lhk>, where 70 = 1 - 7EfcP 2 ( fc ) and < (7! < EfcP 2 (fc)] -1 . The problem 
(|54p now reduces to a quadratic equation for F(k), of which the nonnegative solution is 

F(k) = \y + \^y 2 + A 1P {k)k/{k) (72) 



l-7(p(fc)) / 2fc v 

<fc) \y+y/y*+4rfp(k)k/(k}/ 1 ' 



The right-hand side of (|73p decreases monotonically from [1 — 'j{p(k))]{\/k/p(k)) /y/jffl) a ^ 
y = to zero as y — ► 00, so ([73)1 always has a unique non-negative solution y. One also quite 
easily established the useful bounds 

(which are seen to become tight both for 7^0, where y = 1, and for 7 — > (p(k))~ 1 , where 
y = 0). The matrix (I50p will always have an eigenvalue A = 0, corresponding to the eigenspace 
ek = for all k > 1. The eigenvectors of ([50]) with nonzero eigenvalue A are seen to be 

1 
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where A then follows upon solving 

Upon inserting (|72p this equation takes the following explicit form, with y (which is itself not 
a function of A) to be solved from 



A 



/ Ak(k — \) \ 

1 = l 1 -7(p(*))]( A ^j !/ + ^ y 2 +4:lp ( k ) k /^2 _ 4 1P (k)k(k-l) / 

The right-hand side diverges to — oo for A { and decays to zero for A — > oo. Furthermore, 
provided there exists k > 1 with p(k) > 0, it has singularities for each k with p(k) > at the 
special values A = X c {k), where 

A (k) = ^P(k)Hk-l) , , 

C{ ' (k)[y + y/y*+4ryp{k)k/(k)]* ~ { ' 

Let us define max fc p ( fc ) >0 X c (k) = \ c (k*). We know that the right-hand side of (|77|) decreases 
monotonically on the interval [X c (k*),oo) from oo down to zero. Hence there is always a 
positive solution A of (|77p. and the largest solution X max (Q,p) lies in [A c (fc*), oo). Furthermore, 

A max (Q, P) > max XJk) = max ^ p(k)k(k-l) 

maxW ' J ~ k cK J k (k)[y + y/y 2 +4 1P (k)k/{k)} 2 V ' 

A second simple but effective bound on solutions A > is established easily: 

1-7(pW) / 4fc(fc-l) 

(k) \ [y + y^T^fpi^k/ik)} 2 - 4rf P (k)k(k-l)/X(k) 
l-7(p(fc)) / 4fc(fc-l) v l-7(p(fc)) / Mfc-l) \ 

(A;) \[!/ + ^T%(fc)p)r " W \y 2 +^p(k)k/(k)/ 

and hence 

ww.^) > n-7W)]( ^ra_ ) (80) 

For the trivial choice 7 = 0, where Q(k,k') = 1, we recover the correct results for the 
non-deformed ensemble, viz. F(k) = 1 and lim 7 _»o ^max(Q)P) = {k 2 }/(k} ~ 1- Here 
the second bound (|80p is satisfied with equality. In the opposite limit 7 — ► (p(k))^ 1 , 
where Q(k,k') -> 5 kk >/{p(k)), we obtain F(fc) = p{k)k / ^ (p(k)) (k) + 0(e) and y = 
e(yjk/p(k)) I yj (k) I (p(k)) + C(e 2 ), with e = l—j(p(k)). Our equation for A thereby becomes 

1 = / fc(fc-l)(p(fc)) v 

\p(A;)fc(A-fe+l)/e + X^/pJkjk(p(k))(^/k/p(k)) - p(k)k(X-k+l) + 0{e 2 )' 1 ; 

It follows that for 7 — » (p(k)}^ 1 all nonzero eigenvalues are of the form A = fe*— l+0(l-^y(p(fc))) 
with k* £ {1,2,...} such that p{k k ) > 0. The largest such eigenvalue corresponds to the 
largest k* with p(k*) > 0, so lim 7 _ j ./ J ,( fc \\-i A max ((3,p) = fe*— 1. Thus in the latter limit we 
obtain the transition lines corresponding to a regular random graph with degree k*, which is 
consistent with our earlier observation that for Q(k,k') = S^y / (p(k)) our graphs decompose 
into a collection of disconnected regular graphs, one for each degree k that is allowed by p(k). 
Here we find that the first bound (|79p is satisfied with equality. 



4-3. Phase diagrams 



Our order parameter equations apply in principle to arbitrary choices of the bond distribution 
P(J), the degree distribution p(k), and the ensemble deformation function Q(k,k'). Here we 
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will limit ourselves for brevity to the deformation functions analyzed in the previous section, 
and to the binary bond distribution P(J) = |(1 + n)5{J — Jo) + ^(1 — n)5(J + Jo), with 
n G 



■1,1] and Jo > 0. Equations (|55|56p can now be written as 

r??A max (Q,p) + 1- 



T F / J = 2/ log 



P^SG: r SG /Jo = 2/log 



n\max(Q,p) - 1 

\fKns^{Q,P) + 1 



1 



(82) 



(83) 



The P^F transition is the physical one for n > tanh(/3Jo) and the P— >SG transition is the 
physical one for n < tanh(/3Jo), with a triple point at n = tanh(/?Jo). We will consider only 
two types of degree distributions, both with average connectivity (k) = c: 



Poissonnian : 



power law : 



p(k) 
p(k) 



c /k\, 

c C(3+a) 



$ko + (l-4o) 



ck 



-3-a 



C(2+a) 

Here £(x) denotes the Riemann zeta function £(a: 



C(2+a) 



(k 2 )/c 
(k 2 )/c 



c+l 
C(l+a) 



(84) 
(85) 



J2k>o k 



C(2+a) 
and we take a £ [0,1] to 



ensure that c = (fc) exists (limiting ourselves to c < £(2+a)/£(3+a), so that p(0) > 0, which 
means that c will remain modest), but with the possibility to take the scale-free limit a — > 0. In 
practice, however, in calculating averages over p(k) numerically one has to truncate the values 
of k; here we used k < /c max = 10 8 . For Poissonnian p(k) this has no noticeable implications, 
but for power law p(k) the slow divergence of J2k ^ ~ ^max manifests itself in transition 
temperatures for a — > that should have been infinite but are finite. On the other hand, in 
any finite real system or simulation one will have k < N, so one expects to see also there the 
same effects of bounded degrees (e.g. finite transition temperatures). The 'ideal' situation of 
unbounded degrees and truly scale-free graphs is never realized in practice. The power-law 
distribution ([85]) has the property p(k)\ c = cp{k) c= \ for k > 0. Hence for any function tf){k) 
with i/j(0) = one will have (ip(k)) = c(ip(k)) c= i. As a consequence one finds immediately 
upon checking the various formulae of the previous section that the bifurcation lines for type 
I and type II ensemble deformations are completely independent of the connectivity c for 
power-law distributed degrees. For type III deformations this is not the case. Note, finally, 
that there is no point in choosing regular graphs p(k) = <5j. c , since there the function Q(k, k') 
is always equal to one due to the normalization requirement J2kk' p{k)p(k')Q(k, k') = 1. 

We will compare phase diagrams for the previously analyzed families of deformation 
functions Q(k,k'), viz. the separable ones, the additive ones, and the binary ones. In the 
separable case (type I), where one always has the simple eigenvalue X mSLX (Q,p) = (k 2 }/(k) — 1, 
we have fully explicit expressions for the transition lines that are identical to those describing 
non-deformed ensembles with degree constraints only: 



Poissonnian p(k) 
power law p(k) : 



r F /j = 2/io g 
r F /j = 2/io g 

T SG /Jo = 2/log 



nc+1 



nc — 1 



T SG /Jo = 2/log 
+ (l-Tj)C(2+a) 



y/c-1 



nC{l+a) - (l+7?)C(2+a) 



(86) 
(87) 



VC(l+«) - C(2+a) - y/((2+aj 

(in non-deformed graphs of the type considered here, the transition temperatures for power- 
law distributed degree distributions are independent of the average connectivity). Clearly, in 
non-deformed Poissonnian graphs we can only have an SG phase if c > 1 and an F phase 
if c > 1/ry, whereas in non-deformed power law graphs we can only have an SG phase if 
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c c c 



Figure 1. Continuous bifurcation lines for P^SG (dotted) and P— »F (solid, with 
i] € {0.25,0.5,0.75,1} from bottom to top), for type II deformed ensembles (with 
Q(k, k') = [g(k) + g(k')]/2(g(k))), P( J) = ±(1 + g)5(J - J ) + ±(1 - n)5(J + J ), and 
Poissonnian degree distributions p(k) — c k e~ c /kl. Left to right: g(k) E {1, k, k 2 }. The 
left picture represents the non-deformed ensemble, to serve as a reference. The effect 
of a deformation with Q(k, k') — [k rn + (k') m ]/2(k rn ) in graphs with Poissonnian p(k) 
is seen to be a slight reduction of all critical temperatures with increasing m. 




T_ 

Jo 





g(k) = k 







T_ 

Jo 




a 



a 



a 



Figure 2. Continuous bifurcation lines for P^SG (dotted) and P— >F (solid, with 
g € {0.25,0.5,0.75,1} from bottom to top), for type II deformed ensembles (with 
Q{k, k') = [g(k) + g(k')]/2{g(k)}), P(J) = §(1 + g)S(J - J ) + |(1 - g)S(J + J ), and 
power law degree distributions p(k) ~ fc _3_a . Left to right: g(k) G {1, k, k 2 }. The left 
picture represents the non-deformed ensemble, to serve as a reference. The effect of a 
deformation with Q(k, k') = [k rn + (k') m ]/2(k m ) in graphs with power law p(k) is now 
seen to be a dramatic reduction of all critical temperatures with increasing m. 



C(l+a)/C(2+a) > 2 (giving a < a c w 0.479) and an F phase if C(l+a)/C(2+a) >l + l/n. We 
will not show these lines describing the non-deformed ensembles in a separate figure, but will 
include them as a benchmark when showing data for the type II and type III deformations, 
since in type II models the non-deformed ensemble is recovered for the special choice g(k) = 1 
whereas in the type III models it corresponds to 7 = 0. 

For the additive functions Q(k,k') (type II) the eigenvalue A max (Q,p) depends in a 
nontrivial way on Q(k,k') and p(k), and must be solved from (|70I71|) numerically. Here 
we will choose either g(k) = k or g(k) = k 2 in the function Q(k,k'). For Poissonnian 
distributed degrees this gives G(k) = k/c and G(k) = k 2 /c(c + 1) respectively. For the power 
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c c c 



Figure 3. Continuous bifurcation lines for P^SG (dotted) and P— >F (solid, with 
i] £ {0.25,0.5,0.75,1} from bottom to top), for type III deformed ensembles (with 
Q(k,k') =7o + 7<W)> P(J) = \{l + rj)8(J- J Q ) + \(l - rj)8{J + J ), and Poissonnian 
degree distributions p(k) = c fc e" c /fc!. Left to right: j(p(k)) € {0,0.4,0.8}. The left 
picture represents the non-deformed ensemble, to serve as a reference. The effect of a 
deformation with Q(k, k 1 ) = 70 + 7<5fcfc' in graphs with Poissonnian p(k) is seen to be 
a significant increase of all critical temperatures with increasing 7. 




C C C 



Figure 4. Continuous bifurcation lines for P^SG (dotted) and P^F (solid, with 
77 G {0.25,0.5,0.75,1} from bottom to top), for type III deformed ensembles (with 
Q(k, k') = 70 + 7***'). P ( J ) = + fM- 7 - J o) + |(1 - »7)<y(^ + Jb), and power law 
degree distributions p(k) ~ fc _3_Q with c = 1. Left to right: j(p(k)) G {0,0.4,0.8}. 
The left picture represents the non-deformed ensemble, to serve as a reference. The 
effect of a deformation with Q(k, k') — 70 +7<5fcfc' in graphs with power law p(k) is now 
seen to be a dramatic increase of all critical temperatures with increasing 7. 



law distributed degrees one finds G(k) = k/c and G(k) = k 2 C(2 + a)/c((l + a), respectively. 
Upon solving (|7QI7ip numerically for Poissonnian graphs, we obtain the bifurcation lines as 
shown in Figured! We also show the lines for the non-deformed case g(k) = 1, as a benchmark. 
The deformation causes only minor changes to the phase diagram, mainly a slight reduction 
of all transition temperatures for small values of the connectivity c. When applied to graphs 
with power law degrees, in contrast, the impact of the deformation is much more drastic, as 
shown in Figure [21 This can be understood mathematically on the basis of equations (I70|7ip . 
If we consider the case a — » we only need to inspect what happens to the divergent sums 
over k: one finds for both g(k) = k and g{k) = k 2 that A max (Q,p) ~ ^/log A; max as a | (rather 
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than A max (Q,p) ~ log/c max , as was the case for the non-deformed ensemble). 

Finally, we have solved numerically equations (|73|77p for the case of the binary functions 
Q(k,k') (type III), for r y{p{k)) G {0,0.4,0.8} (with the first value 7 = 0, the non-deformed 
case, serving as a benchmark) and c = 1. Positive values of 7 imply increased connections 
between links with identical degree, which favours especially the formation regular graphs with 
large values of k. Here we always observe a significant increase of all critical temperatures, 
both for Poissonnian and for power-law distributed graphs, see Figure HI The effect becomes 
stronger as c increases. Choosing negative values of 7, i.e. discouraging the formation of links 
between nodes with identical degree, is found to decrease all transition temperatures. Note 
that without the degree cut-off k mSiX = 10 8 , one would have diverging critical temperatures at 
all c > and all a > in the limit j(p(k)) — > 1. 



5. Discussion 



The rationale behind studying interacting particle models on complex random graphs is that 
the latter can be used as solvable proxies for models on specific graphs of a topology for 
which either no exact solution is available (e.g. spin models on cubic lattices), or on which we 
lack precise information (e.g. proteomic networks). We then have to choose an appropriate 
ensemble of random graphs, which is sufficiently simple to allow for analytical progress, while 
incorporating as much as possible of the topology of the specific system one aims to understand. 
Specifying just the degree distribution p(k) of a complex connectivity graph for an interacting 
spin system will clearly not yet permit reliable predictions on the system's phase diagram. 
For instance, the critical temperature of the D-dimensional Ising model on a cubic lattice, 
where p(k) = Sk y 2D, is different from that of a regular random graph with p(k) = 5fc,2iJl- 
The question is then which further topological information on a graph beyond p(k) could be 
added to reduce the entropy of the underlying graph ensemble and make more specific and 
more accurate predictions of phase transitions, while at the same time maintaining the vital 
property that the resulting spin models can be solved analytically. In this paper we have 
established that the proposed deformation of random graph ensembles can be a useful step 
in this direction: it generally allows us to differentiate between models with the same p(k) 
(which can be chosen freely) but different microscopic realizations of these degree statistics, 
the resulting models are still solvable, and its impact on the transition lines can be non- 
nce ligiblfl In practice, when seeking to model a complicated real system with some specific 
given interaction graph c* (and hence a known set of degree {k\, . . . , k^} and a known degree 
distribution p(k)) by a solvable system on a random graph, we could now incorporate at least 
some of the extra topological information by using our ensemble ([2]) with constrained degrees 
ki = k\ for all i, and with a function Q(k, k') that is taylored to the graph c*. This can be 
done by maximizing the log- likelihood of c* for the ensemble ([2]), i.e. by minimizing over Q 
(subject to J2kk' p(k)p(k')Q(k, k') = 1) the quantity 



^[Q] = llogZ JV -l^log 

i<j 



^Q(k*,k*)5^ + (l-^Q(kt,k- 3)r < :s) 



\ For a ferromagnetic Ising model on a square lattice in D = 2 one has Onsager's famous result 
T c /J = 4/log[(\/2 + 1)/(V2 - 1)] w 2.26919, whereas T c /J = 2/log2 « 2.88539 in the degree-4 
regular random graph. Yet both models have the same degree distribution p(k) = #^4. 
§ The only graphs where ensemble deformation is not possible are the regular graphs, with p(k) = 5k c , 
where the constraint Ylkk' P(k)Q(k> k')p(k') = 1 leaves only the trivial choice Q(k, k') = 1. 
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= z + \{k)-\{k) \o g [{k)/N\ - 1 y, 4 lo s k j) + °( Ar_1 ) 

i<j 

= const + log F(fc|Q) - 1 £ c*- log fc*) + ©(AT 1 ) (89) 

where F(fc|Q) is the solution of 

F(fc) = (ky^Pik'WQik, k^F-^k') (90) 
fc' 

This will be the subject of a subsequent study. In addition one would like to study certain 
technical aspects of the present model in more detail, such as the precise physical meaning 
of the function F(k), and the impact of possible replica symmetry breaking (RSB). In the 
present type of model RSB does not change the locations of the P^F or P^SG transition 
lines, but will alter the nature of the solution in the ordered phases and the location of the 
F^SG transition line. 
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Appendix A. Joint distributions of degree and clustering coefficients 



To characterized a graph's local topology we can define for each vertex i the degree fcj(c) = 
J2j c ij (the number of links to this vertex) and the number of length-three loops going through 
this vertex, as measured by rj(c) = J2jk c ij c jk c ki- The clustering coefficient Cj is then given 
by Ci = ri/ki{ki — 1). We write their joint distribution as P(k,r\c) = N^ 1 Ei ^fc,fci(c)<5r,r-i(c)> 
and the asymptotic expectation value of this distribution over the ensemble (|2J) as 



P(k,r) 



lim — — — 

N^oo NZ N 



E E = / ^ r JJm i E KkAW (A.l) 



with 



■>N 



(A.2) 



It will turn out that here we have to expand to higher orders in N than in previous calculations. 
We now re-name all links to/from site i as Sj = £ {0, 1}, while writing all those that do 
not involve site i as Tj£ £ {0, 1}, where j, i G {1, . . . , % — 1, % + 1, . . . , N}. This gives 



sT 



n 



&Q(k i ,k j )s aj , 1 + (i-&Q(k i ,k 3 : 



,. Sj 
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t<j\ej¥=i 



1 



-"TV 



—Q(ki,kj)S SjA + (1-— Q(fci,fej))<y Sj) o 



-21(^8^ + 0^) 



(A.3) 



(A.4) 



At this point we are led to the introduction of the observables 
W sk ((f>) = j^—^Y^8 SSj 8 kk ^{<l>-4>j) 

Clearly E s e{o,i} Efc>o fd<f> W s k(4>) = 1- We also introduce the short-hand p k = N~ 1 J2i$k,k i 
(viz. the empirical degree frequencies, which will only be identical to p(k) for N — > oo). Upon 
introducing the W sk {(f) in the usual manner via suitable 5-functions we can then write 

Pity) = J_ r^e 1 ^ /{dWW} e^-^E.^J^'^^^O^^^O+o^- 1 ) 



x e 



X e -|(fc) 2 E s ' fc ' s » fe »/ d ^'^"^' fe 'W')^» fc »(0")Q 2 (fc',fc")[e" i(2 ^ s '-I] 2 
X e"^ fc ) E s ' fc ' / d4,'W s , k ,W)Q{k>,k>)[ e -^' +*')-!] 



n{E 

1 /"* (ty: 



* d^r(fc) 



-7T 27T 



iV 



Q(ki,kj)S S:1 + (l-^-Q(ki,kj))6 8fi 



-^s+iWkj-W^U'W 



JL r^ e i#, / {d WW} e^-^E, 
Z c 27r 7 



'k' 



x gi(iV-2)(fc> £ s , fcV , fc „ /d0'd0"W s , fe ,(0')W s » fc »(^")Q(fc',fc")[c- i ( 2 ^ s ' s "+^'+^")-l] 
X e -3 (fc)2 E s ' fc V' fe » / d0'^"W s , fe ,(0')^ fc »(0")Q 2 (fc',fc")[e- i(2 "' s ' s " + ^ + ^" ) -l] 2 

x e" 



i <*> E s ' fc ' / d^Wi/ fc #(^)Q(fc',fc')[e- 



»n [■ 

3 + 



(i 7^ Q(ki,kj) 



J —IT 



2tt 



(A;) f' r di'e l[ ** r ' Vu i ( * H 



J_ r WfMki [{dWdW} e i ^- 1 )E sfc / d ^^fe^')H4 fc (^')-<fe>E fc PfeQ( fc - fc )+o(JV- 1 ) 
i(JV-2)(fc> £ s , fcV , fc „ / d^d0"w s , fc ,(0')^» fc „(^")Q(fc',fc")[c- i ( 2 ^ s ' s " + ^'+^")-i] 



-'C 

x e2 
x e~3 



\(k) 2 E s ' fc ' s '' fe » / d0'^"W s , fe ,(0')^ fc »(0")Q 2 (fc',fc")[e- i(2 " ,s ' s " + ^ + ^" ) -i] 2 



xe "K fc )E s ' fc ' /d^M;, fc ,(^)Q(fc',fc')[e- 2i W''+*')-l]+ArE fc Pfcl«8/!^e i[ *'*- ,i ' 0fc( *' )] 

exp (fe) V Pfc Q(fc i} fc)^^ s log / ^fiil^fci-^^)] 

P I N '4^*^ l ' ' J^d^e^-^W')] 



2vr 



= _}_ T # e i0fci /"{dM^dlU} e Ar *( w '^)+*W^^) +n (^' A; "^ +c '( Ar_1 ) (A.5) 

7-7T 27T 7 



with 



tt(w;w^) = i]T Td0 w ak (tf,)w sk (tf>) + 5> fc io g f ^e^-^wi 
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+ \( k ) E /#d<m*(»*^OQ^ (A.6) 

sks'k' 

$(W,W,if>)= — (k) J2 y'd0d0 , W afc (0)^ fc /(0 , )Q(fc.fc')[e" i(a0 *'' + * + *' ) -l] 

sks'k' 

~ \( k ) 2 E / ^W 8k (<f>)WMW^ k ') [e -^ss' +< p + <p') _ 1]2 (A _ 7) 

sks'k' 

sk sk 

r f 71 " d^V^* - ^ 1 *^) - ^ i r w d6' 
Q(W,ki,6) = (k)Y 'p k Q(ki,k) 1 - log / — 



(Ai 



Using the normalization identity Pj(0) = 1 we may then also write 

P«M) = J- » » ; (A.9) 

J^d^ e 1 ^ /{dWdiy} eN^{w,Wfi)+HW,Wfi)+n(w,k t ,<t>)+o{N- 1) v ; 

and, upon defining P(r|/c) = P(k,r)/p k : 

P(r\k)= £^tj^ LkW (A10) 

ffdWdWT- e iV *( W '^^) + *( w ''^) T d<£ J**+n(ifr.*.0 

L fc (^) = lim — ^ = = = (A.ll) 

KKY > iv^oo jjd^dW} e Nm(w,w$)+®{w,w$) f^ dcj) e i<t>k+n(w,k,<f>) v ; 

We next need to find the saddle-point (s) of the function (|A.6p . by variation of {W, W}. 
Functional differentiation with respect to W and W gives the following equations, respectively: 

iW sk {4>) = - (k)J2 /#' W s/k >(</ ) , )Q(k,k')[e- i ^ ss '+'t>+^-l} (A.12) 

s'k' J 

e i[4>k-W 0k (<t>)] 

W sk {$) = S s0Pk - (A.13) 
f_d<j>' e l \<p k -w 0k {<t>')\ 

Upon eliminating W and defining W sk ((p) = S soPkXk(4>) > we obtain an equation for Xk{4) only: 

e i*fc+(fc)Efc/PV < 3( fc » fc ') e " i * JV xvW-'K^' 

Xfc ^^ ~ ^ e i^fc+<fc>£ fc /Pfc'Q(fc.fc')e- i *'/d^' Xfe 'W")e-^" ^ A ' 14 ' ) 

One defines a k = J d<f>' Xfc(^') e-l< ^ an d bfc = (&) J2k' Pk'Q(k, k')a k >, and finds after some simple 
manipulations that ao = and a^x) = This leaves a closed equation for the b k , which 

shows that limAr^oo b k = (k)F(k), see and a corresponding formula for Xfc(</0 : 

6 fc =(k)J2Pk'Q(k,k')k'/b k/ (A.15) 

fc'>0 

Xfc(0) = ^-fee-^)- fc exp[6 fc e-^ (A.16) 

At the relevant saddle-point, we find as a direct consequence of the form W sk (4>) = $soPkXk{4>) 
that the functions *$>(W,W and <&(W,W,iJj) are both independent of the variable ip. This 
ensures that expression (IA.11I) is well-defined, but it also gives us L k (ip) = 1, and hence 

P(r\k)=5 r0 (A.17) 

We conclude that in our ensemble ([2]) the fraction of nodes in a loop of length three vanishes in 
the limit N — * oo, independent of the degree distribution p(k) and independent of the choice 
made for the deformation function Q(k, k'). 



